Some remarks about the domain statistics in the Ising model 
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In this paper we describe a relation between the Zipf law and statistical distributions for the 
Fortuin-Kasteleyn clusters in the Ising model. It has been shown,that histograms for fixed domain 
masses present the right-skewed distributions. 
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In our paper we consider power law distributions, espe- 
cially the Pareto distributions. These distributions have 
been investigated for the last few decades in connection 
with various subjects for examples languages analy- 
sis of populations in cities|2j,the turnover of Europes's 
largest companies Q, biology [j] and physics [E[ , [1] , [3] , [1] , 

0- 

Although in this paper we consider a problem rather 
well understood qualitatively, it is still far from being 
closed and it requires some explanation. We study the 
size distributions of the Fortuin-Kasteleyn( F-K)clusters 
obtained by the Monte Carlo simulations in the two- 
dimensional Ising model. It is well known, these clusters 
obey the power law at criticality. In our case the 
measure of size is the number of spins in the cluster 
( the mass of the cluster), but it can be also another 
measure of size as a disk covering cluster or an area 
enclosed by boundaries @|. The corresponding fractal 
exponent may be related to the criticality exponent as is 
shown in[lo|.or more recently in [llj . Similar situation 
we have for a percolation. This idea for the Ising model 
in 2d may be found in[T2| . 

The idea of the power law is explaned as consequence of 
the scale invariance at a criticality due to the divergence 
of the correlation length . Away from the criticality, the 
range of lengths, where the system is fractal, is limited 
by the finite correlation lenght, so that the power-law 
tail of the distribution is suppressed. 

We consider the Hamiltonian for the Ising model 



H 



QZ QZ 



(1) 



with a sum over all neighbouring pairs (z-th compo- 
nents) of spins. Usually it is assumed that the crystal 
lattice of a ferromagnet is regular and in each site of a 
lattice the spin is localized with the value S z = 1 or 
S z = — 1. Further we have: 



Jij — 



J if i,j are neighbouring pairs of spins 
if i,j are opposite each other. 



(2) 
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Two spins i and j interact with each other by an 
energy ± J, Sf with — J if both spins are parallel and 
+ J if they are opposite each other. The energy needed 
for flipping of one spin is 2 J. 

Results presented hereafter are obtained by applying 
the Monte Carlo techniques, based on the Swendsen- 
Wang cluster algorithm [13j to the two-dimensional Ising 
model with periodic boundary conditions. In this algo- 
rithm clusters of spins are created by introducing bonds 
between neighboring spins of the same orientation, with 
the probability 1 — exp(— rM|;), where fcs is the Boltz- 
mann constant and AE is the energy required to trans- 
form a pair of equal spins to a pair of opposite spins. To 
each cluster we assign at random —1 or +1 orientation. 
The different clusters are getting independent orienta- 
tions. Starting with the simulation having a random dis- 
tribution of a half of the spins up and a half of the spins 
down and using the Swendsen-Wang algorithm with a 
low temperature one sees the growing domains, in which 
spins are parallel. We have two kinds of domains: with 
spin up and with spins down. 

The phase transition appears at the critical point 
T = T c (T c = 2J/k B ln(l + V2)). The difference be- 
tween the number of spins up and down M is propor- 
tional to the magnetization and near the critical point 
vanishes as (T — T c )&, where for the dimension d = 2, 
P = i. The correlation length £ ~ \T-T C \", while the 
magnetization is proportional to £ _ ^/ ly '. In a finite sys- 
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tern at the critical temperature T c one can replace £ by 
L, hence M ~ L d ~P/ v — L D with the fractal dimension 
D = d-p/v (d<4). 

In this paper we want to show that the right-skewed dis- 
tributions of the clusters are present in the histograms 
with the fixed masses and we explain the relation be- 
tween the phase transition , the Zipf law ,the power law 
and the Pareto distributions. 

Firstly in this paper we want to underline that the Zipf 
law denotes greater restriction than the existence of the 
power law expansion. The power law is presented by the 
Mandelbrot law in the form 1141: 



(3) 



where C and a and p are some constants and x denotes 
the value of the object, k is the rank order of the ob- 
ject. The object with the largest variable value is ranked 
as the first, the next with a smaller value is ranked as the 
second and so on. In this way we obtain the rank order 
k of the object. The Mandelbrot law undergoes the Zipf 
law when a = and p = 1 : 



C/k 



(4) 



This equation performs a straight line in a double log- 
arythmic plot with the slope equal to 7r/4. 
In the case of the Ising model x denotes the domain mass 
of the rank fc(see Fig. [T|). 





1/T=0.44068 




1/T=0.4 




1/T=0.35 




1/T=0.25 




ln(k) 



FIG. 1. The log-log distribution of the domain masses Xk ver- 
sus the rank order index k for the two-dimensional Ising model 
for T = 2.2692(the critical temperature) (0),T = 2.5000 (*), 
T = 2.8571 (□) and T = 4.0000 (o) for L = 1000. 

The straight lines with the slopes near —1/fJi are given 



in Table U correspond to temperatures greater than the 
critical temperature T c (T c = 2.2692). 



TABLE I. The estimation of indexes p with standard devia- 
tions s. 



T 


/" 


s 


4.0000 


6.7136 


0.0477 


2.8571 


4.7018 


0.0603 


2.5000 


3.8240 


0.0695 


T c ~ 2.2692, 


0.9854 


0.0874 



The estimations of the indexes p for the two- 
dimensional Ising model with different temperatures 
based on 1000 realisations for L = 1000 lattice are pre- 
sented in Table |U and Fig.Q] 

It is seen from Fig. Q] that a slope arctan(l/p) of the 
line only for the critical temperature is near 7r/4 for the 
critical temperature and the Zipf law is satisfied. For 
the temperatures above critical temperature, slopes are 
smaller than 7r/4 and the Zipf law is not valid. One can 
notice that the validity of the Zipf law when the phase 
transition is present is in agreement with the relation 
1/p = \2H — 1| [TU, where H is the Hurst exponent.lt 
is well known that the system exhibits infinite long cor- 
relations at the critical point, which corresponds to the 
biggest value of the Hurst exponent. In such a case one 
obtains p = 1. 

The Pareto distribution implies the Zipf law [lTj , [1] ,but 
not inversally - for instance the Zipf law is satysfied also 
in the case of the hyperbolic distributions fhsj . 
The probability density of the Pareto distributions is de- 
fined as follows: 



p(x) = (p/x )(x a /x 



(5) 



for x > xq, where x$ denotes a typical scale. 

One can define so called fi- variable [l],[l6j]. The 
main property of [/invariable is that all its moments 
Tfiq —< X q > with q > fi are infinite. The question 
in classic probabilistic theory answered by Levy and 
Khintchine shows how to characterize the limit distribu- 
tion of the sum of N independent random variables x. 
Suppose that p(x) decreases for a large x such as 
then considering power-law distributions we can dis- 
tinquish the following ranges for the critical exponents p : 

(i)For p > 2, the moment of the first order mi and the 
moment of the second order m-i are finite. 
For this range of /i,the two parameters mi and vni are fi- 
nite and for large N the central limit theorem applies and 
the statistical distribution of the system tends towards a 
Gaussian distribution. 

( ii)For 1 < p < 2 only the moment of the first order is 
finite but the moment of the second order is infinite, 
(iii) For p < 1 the moments of the first and the second 
order are infinite. For a large number of elements of the 
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system it is described according to the Levy stable dis- 
tribution, for instance the Pareto distribution. 
But magnetic domains described by the clusters of thi 
Ising or the Potts models are not independent variabes 
Although we do not have a correspondence between / 
and the slope of distributions described in (i)-(iii), re 
spectively, the presence of the Pareto-type dystrybution; 
are also possible in the case of the dependent variables . 



figurations is 10000. Histograms are presented for masses 
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FIG. 2. The dependence between the number of domains wit] 
mass x normalizedby the number of configurations (N x /N 
from the mass (:r)in the log-log scale, for the two dimensiona 
Ising model near the criticalpoint T c , for 10000 configuration 
and L = 600. 
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FIG. 3. Histogram of the fixed domain mass m = 5(a),m = 
10(b), m = 50(c) and m = 70(d) for the inverse temperature 
j9 = 025, for the system with L = 600, for 10000 configura- 
tions. 



x = 5, x = 10, x = 50 and x = 70 for ,for a system with 
L = 600 near critical point (T c = 2.2692). 
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The straight line of Fig. [5] presents the Pareto distri 
bution. 

In our case, away from the criticality, the Gaussian 
type distributions are not present when /i > 2. In thi; 
region the range of lengths, where the system is fractal, i 
is limited by the finite correlation length and the power 
law tail in the distribution is suppressed. 

Fig. [3] presents the distributions of the number of 
configurations with defined mass. The total number of 
configurations is 10000. Histograms are presented for 
masses x = 5, x = 10, x = 50 anda; = 70 for , for 
a system with L — 600, beyond the phase transition 
inverse temperature /3 = 0.25. 

The test for the Gaussian distribution was performed 
using the Kolmogorow-Smirnov method and the chi- 
square test. The Gaussian distributions for the num- 
ber of clusters of a fixed size are present only for the 
masses smaller than 5. For example for the temperature 
T = 4.000 and the cluster of size 5 the p-value of the chi- 
squared test was equal toO. 0771. For the clusters of larger 
size we obtain the right-skewed distributions. 

Fig. 2] presents the distribution of the number of con- 
figurations with a defined mass. The total number of con- 
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FIG. 4. Histogram of the fixed mass m = 5(a), m — 10(b), 
m = 50(c) and m = 70(d) for the critical temperature, for 
the system with L = 600, for 10000 configurations. 

In this case the right-skewed distributions are 
present. The Gaussian distributions we observe only for 
masses smaller than 3. For example for the cluster of size 
5 the p-value of the chi-square test was equal to 0.01174. 

The main conclusion from this paper is following: 
In the critical region the Zipf law is satisfied and the 
power law distributions, the Pareto distributions are 
present. The histograms of the numbers od domain con- 
figurations, which possess the same number of domains 
present the right-skewed distributions mainly, for the 
critical point as well beyond it. The Gaussian- type dis- 
tributions appear only for small masses . 
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